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ABSTRACT 

We have used realistic Comptonization models to fit high quality BeppoSAX data of 6 Seyfert galaxies. 
Our main effort was to adopt a Comptonization model taking into account the anisotropy of the soft pho- 
ton field. The most important consequence is a reduction of the first scattering order, which produces a 
break (the so-called anisotropy break) in the outgoing spectra. Thus anisotropic Comptonization models 
yield spectra with convex curvature. The physical parameters of the hot corona (i.e. the temperature 
and optical depth) obtained fitting this class of models to broad band X-ray spectra are substantially 
different from those derived fitting the same data with the power law -I- cut-off model commonly used 
in the literature. In particular, our best fits with Comptonization models in slab geometry give a tem- 
perature generally much larger and an optical depth much smaller than derived from the power law -|- 
cut-off fits, using standard Comptonization formulae. The estimate of the reflection normalization is 
also larger with the slab corona model. For most objects of our sample, both models give Compton 
parameter values larger than expected in a slab corona geometry, suggesting a more "photon starved" 
X-ray source configuration. Finally, the two models provide different trends and correlation between the 
physical parameters: for instance, with the power law -I- cut-off fits, we obtain a correlation between 
the reflection normalization and the corona temperature whereas we find an anti-correlation between 
these parameters with the slab geometry. These differences have major consequences for the physical 
interpretation of the data. In the framework of reprocessing models, the cut-off power law best fit results 
suggest the thermal corona to be dominated by electron-positron pairs. On the contrary, the slab corona 
model is in better agreement with a low pair density solution. 

Subject headings: galaxies: active; galaxies: Seyfert; X-rays: galaxies 



1. INTRODUCTION 

The broad band X-ray spectra (2-300 keV) of Seyfert 
1 galaxies are generally well fitted by a cut-off power law 
continuum with superimposed secondary components like 
a neutral iron line and a reflection hump. The latter com- 
ponents have been detected in the end of th e eigh ties with 



al. 



1990 



Matsuoka et al. 



1990; Pounds 



GING A (Pir o et 

et al. |1990D and are considered as signatures of repro- 
cessing of the primary X-ray emission in surrounding cold 
matter. A high energy cut-off in the X-ray/7-ray part of 
the spectrum was first unambiguously detec ted b y SIGMA 
and OSSE in NGC 4151 (Jourdain et al., |1992| ; Maisack 



et al. 1993). A high energy cut-off was also requested to 



fit the average X-ray/7-ray spectra (obtained from data of 
differ ent satellites) of a sam ple of Seyfert (Zdziarski et al. 

195 ; Gondek et al. 1996| ). The broad band capabilities 
of the Italian-Dutch satellite BeppoSAX allowed to detect 
high energy cut-offs in other Seyfert 1 galaxies, about 
half of the fifteen observed with this satellite. Values of 
the hi gh en ergy cut-offs range from 70 keV to ~ 300 keV 



(Matt 2000, hereafter MOO). 



A cut-off power law plus reflection model (the so-called 
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PEXRAV model in XSPEC, Magdziarz & Zdziarski |1995D 
depends on three parameters: the spectral index F, the 
high energy cut-off Ec and the reflection normalization i?, 
the formers characterizing the shape of the primary con- 
tinuum. The physical interpretation of these parameters is 
generally done in the framework of thermal Comptoniza- 
tion mechanism, which is commonly believed to be the 
origin of the X-ray emission of Seyfert galaxies. Using 
approximate relations (cf . section ^) , it is then possible to 
derive, from values of F and i?c, values of the temperature 
kTe and optical depth r of the "comptonizing" hot plasma 
(the so-called corona). 



If such approximations are sufficient in isotropic ge- 
ometries (e.g. a central soft photon source surrounded 
by a spherical hot cloud), where Comptonization process 
produces roughly power law spectra with a high energy 
cut-off, strong discrepancies may appear in anisotropic 
ones, especially for small optical depth and large tem- 
perature. Indeed, the (possible) anisotropy of the source 
of seed photons (as seen by the corona), may introduce 
anisotropics and modifications of the outgoing spectrum. 
The largest effect occurs when soft photons are emitted 



2 



by a plane (disk) below the corona. In this case, photons 
backscattered towards the disk in the first scattering are 
necessarily produced in " head-on collisions" and therefore 
have an energy gain larger than average, while photons 
scattered towards the corona (i.e. in the forward direc- 
tion) have an energy gain smaller than average. As a 
consequence, the contribution of the first scattering order 
to the outgoing flux is significantly reduced producing a 
spectral break (the so-called anisotropy break) in the spec- 
trum. Below the high energy cut-off the spectra is then 
better described by broken power laws than by a simple 
power law c ontrary to what is generally believed (Petrucci 



et al., |2000| , hereafter POO). These anisotropic effects be- 
come however less important for larger optical depth (i.e. 
r > 1) since they are "smoothed" by the larger number of 
Compton scatterings suffered by the soft photons. 

It is also worth noting that the energy of the anisotropy 
break which is, in practice, clo se to t he peak energy of the 
third scattering order (Haardt 1993| ), depends on kT^ but 
also on the temperature of the seed soft photon plasma 
kThh. Consequently, if in the case of isotropic geome- 
try the spectral shape is relatively independent on kTbt, 
which fixes only the low energy boundary of the emitted 
spectrum, in anisotropic geometry the shape of the Comp- 
tonization spectra depends on it. 

In a recent paper (POO), we have applied accurate Comp- 
tonization models to fit high quality data of the BeppoSAX 
long look at NGC 5548. We have underlined the impor- 
tance of the anisotropic effects and, in particular, we 
have shown that due to the presence of the anisotropy 
break, above which the intrinsic continuum is steeper, the 
temperature of the hot electrons estimated from Comp- 
tonization models in anisotropic geometry can be much 
larger that the one derived fitting the data with a power- 
law -I- cut off model. 

Besides, we have shown that these different models 
may provide different (and even opposite) relationships 
between physical parameters. For instance, still in the 
case of NGC 5548, the spectral softening of the spectrum, 
observed during a small flare in the central part of the ob- 
servation, goes with an increase of the high energy cut-off 
(i.e. of the corona temperature), in the case of PEXRAV, 
but with a decrease of the corona temperature when using 
an anisotropic Comptonization model. 

In the same spirit of POO, we extend our application of 
accurate Comptonization models to a larger sample of ob- 
jects. We report in this paper the results of this study. 
The paper is organized as follows. In section |2| we briefly 
present the main results of the BeppoSAX observations of 
Seyfert galaxies. We also present the sample of objects we 
use for the purpose of this paper. We describe the mod- 
els we used and the fitting procedure in section |3l We 
also detail in this section the fit of the soft excess/ warm 
absorber (hereafter WA) features observed in each object. 
We present the results of the fits in section ^ and discuss 
their physical interpretations in section ^. We then con- 
clude in the last section. 

2. THE DATA 



2.1. Overview of the BeppoSAX observations of Seyfert I 

galaxies 

Two major BeppoSAX Core programs were devoted to 
classical Seyfert 1 galaxies aiming to study the broad band 
spectra (PI: G.C. Perola) and spectral variability (PI: L. 
Piro). Thirteen sources have been observed on aggregate 
so far, since the launch of the satellite in 1996. 

Reviews summarizing the main results on t hese sources 
have been presented by Piro et al., ( 1999D and MOO. 
For all objects, the spectral models used by the au- 
thors included a power law with exponential cut-off (i.e. 
fs oc E^^^ e^-^^-^") characterized by a photon index P 
and a e-folding energy Ec, a gaussian to fit the neutral 
Iron line near 6.4 keV and a Compton reflection contin- 
uum characterized by a reflection normalization R. Soft 
excess and/or WA features were added if required by the 
data. 

The major findings are the foUowings: 

• A high energy cut-off is generally observed, with an 
e-folding energy typically in the range 70-300 keV. 

• All sources show the presence of a reflection compo- 
nent, so confirming the G INGA and ASCA results 
(Nandra & Pounds, |1994|). 



• Only two objects (Mkn 509 and Mkn 841) show the 
presence of a soft excess. This sharply contrasts 
with the result s obtained with EXOSAT ( Turner 
& Pounds |1989| ) or ROSAT (Walter & Fink, |1993D 
where most of the Seyfert 1 observed seemed to show 
suc h exce ss (for a discussion of this point see Piro et 



al., 1997). On the other hand, WA features (mostly 



the O VII and O VIII edges) are generally observed 
in the LEGS range. 

• A correlation between the power law index P and 
the amount of reflection component R is indicated 



in the same sense as found by Zdziarski et al. (1999 



hereafter Z99). However, the issue of the reality of 
this correlation is still open since it may be due to 
the fact that the two parameters are strongly related 
in the fit ting p rocedure (see for example Vaughan & 
EdelsonpOOq). 



2.2. Our sample 

From the sample of Seyfert 1 galaxies observed by 
BeppoSAX , we have selected observations with good 
signal-to-noise ratios at high energies as measured with 
the PDS instrument. We therefore favored objects with 
hard X-ray spectra. Such observations are best suited for 
an analysis of the high energy spectra. 

The subsample of objects studied in this paper is re- 
ported in Table 0, with the name of the different objects, 
the date of the BeppoSAX observations (some of these ob- 
jects have been observed several times), the exposure time 
and the statistical significance in the PDS instrument. 
We also present results from the observation of ESO 141- 
G55, even if its detection in the PDS instrument is rel- 
atively weak, since this source was recently observed for 
the first time by BeppoSAX . 
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NGC 4151 was observed three times by BeppoSAX , 
two times in 1996 during the SV-phase (in July and in 
December), and the last one in January 1999 during A02. 
However during the observation of July 96, the source 
experienced a strong increase of the 2-10 keV flux by a 
factor two on time scales of one day, wi th sign ificant spec- 
tral variability as shown by Piro et al. ( 2000a ). According 
to these authors, we choose to study separately the low 
and high state occuring during this period (and indicated 
as 96L and 96H in the following). Unfortunately, during 
the high state, the LEGS was off. Only the MEGS and 
PDS data are therefore available for this state. We have 
not used, in this paper, the observation of December 96 
since it does not differ substantially from the low state of 
July 96. 

Goncerning Mkn 509, its observation has been split in 
two parts during the A02. Since no significant spectral 
variability was observed between the two data sets (Per- 
ola et al. |2000| ), we will only refer here to the combined 
spectrum. 

In this paper, we will be concerned with data from 
three of the four instruments on board BeppoSAX : the 
Low Energ y Goncentrator Spectrometer, LEGS (Parmar 
et al. 1997) covering the 0.15-10 keV range, the Medium 
Energy Goncentrator Spectrometer, MEGS (Boella et al. 
1997) covering the 2-10 keV range an d the Phoswich de- 
tector system, PDS (Frontera et al. 1997) covering the 
range 13 - 200 keV. Due to some problems that remain in 
the spectral analysis of LEGS data above 4 keV, we shall 
use only data in the range 0.1-4 keV. 

LEGS and MEGS event files and PDS pha files were 
downloaded from the BeppoSAX SDG archives. The spec- 
tral counts were extracted from a circular region of 4' and 
8' radius around the source centroid in the MEGS and 
LEGS images respectively. We used the data of the three 
(or two, for observations done after May 1997) MEGS 
units merged together to increase the signal-to-noise ra- 
tio. In XSPEG (Arnaud 1996), we use the last updated 
responses (Sep. 1997) and background matrices of each 
instrument. 



3. MODEL FITTING 

3.1. The primary continuum 

We fit the data using two different models for the 
primary continuum: 1) an exponentially cut-off power 
law plus a reflection component from neutral material 
(PEX RAV model of XSPEG, Magdziarzk & Zdziardski 
1995) and 2) a thermal Gomptonization spectrum from a 
disk-f corona configuration in slab geome try. The latter 
was obtained using the code of Haardt (1994, hereafter 
we will refer to this Anisotropic Gomptonization Gode 
as AG2). This code derived the angle-dependent spectra 
of the disk-corona system using an iterative scattering 
method, where the scattering anisotropy was taken into 
account only in the first scattering order. It includes also 
a refiection compo nent described following White, Light- 
man & Zdziarski ( |198^ ) and Lightman & White ( |1988D 
and assuming neutral matter. The spectral shape of the 



reflected photons is averaged over angles. It is multiplied 
by a first normalization factor which depends o n the in- 
clination angle (see Ghisellini, Haardt & Matt |1994| for 



details). In addition, the usual R normalization is left free 
to vary in the fit procedure, so that, for the given inclina- 
tion angle, R = 1 corresponds to a solid angle, subtended 
by the reflector, of 2tt . We do not apply Gomptonization 
models for other geometries (e.g. hemisphere) since, as 
shown for NGG 5548 (POO), the results follow closely those 
obtained with the slab geometry. 

The fit parameters of AG2 are the temperature of the 
corona fcTe, its optical depth r, the temperature of the disk 
/cTfib (assuming a black body soft emission) and the reflec- 
tion normalization R. On the other hand, the PEXRAV 
continuum depends only on 3 parameters: the e-folding 
energy of the cut-off power law Ec, the photon index T 
and the reflection normalization R. 



3.2. Soft X-ray excess and Warm Absorber 

Most of the objects present WA features and/or a Soft 
excess in the low part (i.e. in the LEGS energy band) of 
their spectra. Good and realistic fits require to take these 
features into account. Detailed analyses of these different 
components are however beyond the scope of this paper, 
which is instead focused on the high energy continuum 
of the sources. Then, as a first approximation, we have 
added simple components (edges, gaussian, blackbody...) 
to the primary continuum to reproduce the main features 
present in the soft part of the spectra. We discuss briefly 
the case of each source in the following. 



3.2.1. NGC 4151 

The BeppoSAX data of this galaxy have already been 
analyzed by Piro et al. (2000a, 2000b), assuming a simple 
cut-off power law for the primary underlying continuum. 
A promine nt soft excess, known to exist in this object 
(Holt et al. |1980D below the --2 keV cutoff due to a high 
absorbing column, is clearly detected in the observations 
of 1996 (low state) and 1999. It is generally modeled by 
allowing some fraction of the central source to be com- 
pletely uncovered or viewed through a lower column (the 
so-called "leaky absorber"). In addition, there is evidence 
of a possibly separate soft excess component that does not 
vary with the 2-10 keV continuum (Pounds et al., 1986; 
Perola et al., 1986). Part of this emission comes from 
an extended X-ra.y region know n to e xist in this object 
(Elvis et al. 



1983 



Morse et al. 1995). Recent Ghandra 
observations with the High Energy Transmission Grating 
Spectrometer have revealed detail spectra of this source 
(Ogle et al., 2000[ ) and allow for the first time precise 
measurements of the ionization, temperature and kinetics 
of the extended soft X-ray emission region. 

For the scope of the present paper, we model the com- 
plex soft X-ray emission i n a (relatively ) simpler manner. 
According to Piro et al. (2000a 2000b), we have used a 
dual absorber that is the source is covered completely by 
a medium with a column density iV^fi and with a covering 
fraction fcov by a second medium with a column density 
Nhcov On the other hand, the soft excess below 2 keV 
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was modeled by a scattering component plus an ultrasoft 
component described by a thermal bremsstrahlung, both 
absorbed by the galactic column density. 

Due to the lack of LEGS data for the high state of 1996, 
we cannot constrain the different components of the soft 
X-ray range. We thus analyze the data only above 3.5 keV 
adding a simple photo-electric absorption characterized by 
a column density Nfji. 

We have reported for the three observations our best fit 
parameters Nhi, Nhcov and fcov, obtained with PEXRAV 
to model the continuum, in Table ||. We note a significant 
increase of N hi between 96 L and 99. As already noted 
by Piro et al ( |2000a , 2000b), such variations suggest that 
the structure of the dual absorber has slightly changed 
between the two observations. 



3.2.2. NGC 3783 

Evidences for deep absorption features in the 0.1-1.5 
keV band of NGC 3783 have already been detected in 
ROSAT (Turn er et al. 



1993) and ASCA observations 



(George et al. 1995), the main ones being the O VII and 
O VIII edges at 0.74 and 0.87 keV respectively. More re- 
cently, in January 2000, NGG 3783 was observed using the 
High Energy Transmission Grating Spe ctrom eter on the 
Ghandra X-ray observatory (Kaspi et al. , 2000| ). These au- 
thors detected a large number of emission and absorption 
lines, strongly confirming the presence of warm absorb- 
ing/emitting media embedding the central X-ray source. 

An emission feature near 0.5 keV has also been claimed 



in ASCA spectra of this source (George et al. 1998) and 
has been interpreted as an O V II emission line (expected 
at 568 eV). Kaspi et al. ( |200C| ) clearly see this emission 
feature in the Chandra spectra, confirming for the first 
time the detection of this component. 

According to these results, we have included in our 
models two absorption edges, near 0.74 and 0.87 keV, and 
a gaussian peaking near 0.5 keV to fit the O VII line. The 
best fit results are reported in Table ||. As suggested by 
the recent observation of Chandra, the edge near 1.3-1.4 
likely corresponds to contributions by Fe L and Ne K 
edges whereas the second one, centered near 0.8-0.9 keV, 
may be due to the O VIII edge. 

Concerning the line, we have checked that its addition 
significantly improves the fit with Ax^ = 16 (according to 
the F-test, this difference implies a probability less than 
0.01 that the line is due to random fluctuations in the 
counts). This feature is thus detected by BeppoSAX with 
a high confidence level. The corresponding best fit central 
energy of the line was Eune —0.55 keV with an EW of 
~100 eV. Such EW is large, and it is unlikely that this 
emission features was only produced by O VII emission 
lines. In fact, in the more extreme case where we see the 
blends of the three O VII lines, we expect an EW of only 
30 eV (De Rosa, private communication). The large EW 
and small energy of the line we obtained can be due in part 
to a bad modeling of the continuum near this component. 



3.2.3. Mkn 509 

Mkn 509 was observed by BeppoSAX in 1998. The 
observation was split in two parts during the A02, the 
first one being done in May 1998 and the second one in 
October 1998. As already said, a hardness ratio analysis 
of the counts shows marginal evidence of spectral varia- 
tions between the two observations. They have thus been 
combined and we applied the spectral analysis to the in- 
tegrated counts. 

These da ta ha. ve already been studied in detail by Per- 
ola et al. (|2000[ hereafter PEROO). Assuming a cut-off 
power law model for the intrinsic continuum of this source, 
the data show at low energy the presence of additional, 
superimposed features due to O VII and/or O VIII edges 
as well as a strong soft excess. PEROO described the latter 
by a "soft" power law component {AgE"^''). They then 
obtained a significant improvement of the goodness of the 
fit (Ax^=20), the best fit values for As and Tg being 
1.11x10"^ cm~^ sec~^ kev"-'^ and 2.5 respectively. 

PEROO underlined the fact that the addition of this 
"soft" component had also the interesting effect of pro- 
viding a more astrophysically sound global solution, with 
a reflection normalization R 0.6 and Iron line width 
(jpe ~ 0.4 keV (compared to R ~2 and ape ~3 keV with- 
out this soft component), giving some confidence on the 
correctness of the model. 

Following these authors, we have added this "soft" 
power law to our fit, fixing A^ and F^ to the best fit values 
given above. We also add an edge with the best fit energy 
and optical depth found by PEROO, i.e. Eedge = 0.74 keV, 

Tedge = 0.06. 

3.2.4. ESQ 141-G55 

The presence of a soft excess was already sugges ted fo r 
this source from ROSAT observations (Turner et al 1993). 
However the BeppoSAX data do not show any evidence of 
such component nor WA features. It could be partly due 
to the low S /N of this observation. Consequently, our fits 
were done without any additional components in the soft 
band other than the absorption by the galactic column 
density (cf. Tables | and |) . 



3.2.5. NGC 5548 

The WA observed in this sourc e was already discussed 
in POO and Nicastro et al. ( 2000| ) and the reader can re- 
fer to these papers for more details. For the purpose of 
this work, we adopt the results of POO who have added 
two edges, at Ei^ 0.74 keV and £;2=0.87 keV, to fit the 
Oxygen absorption features observed in the LEGS band. 
The corresponding optical depths were fixed to the best fit 
values obtained with PEXRAV, i.e. ri=0.51 and T2=0.12 
respectively. 

3.2.6. IG 4329A 

The X-ray spectrum of IG 4329A is known to display 
strong neutral absorption, below 3 keV, due to a gas 
column about ten times greater than the g alactic value 
Nh.^ 



4.5 X 10^0 cm-2 (Elvis et al. 1989). This gas is 
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supposed to be associated with the disk of the host 
which is oriented nearly edge-on (Petre et al. 1984) 



alaxy, 



At low energies (below 1-2 keV), ROSAT and ASCA ob- 
servations have also shown evidence of two rather strong 
edges, one consis tent with O VII, t he oth er with O VIII 
(Madejski et al. |1995| ; Cappi et al., |1996| ; Reynold |1997 



George et al. |1998| ) 



The BeppoSAX observation of I C 43 29A has been al- 
ready analyzed by Perola et al. (1999) using PEXRAV 
as primary continuum. To model the soft range of the 
spectrum, they added a neutral column of gas in addi- 
tion to the galactic one and included two absorption edges. 
This model gives a good fit to the data and we adopt the 
same description in our fits. We fix the energy and opti- 
cal depth of the two edges and the neutral column densit y 
to the best fit values obtained by Perola et al. ( |1999| ), 
i.e. El =0.73 keV, ti=0.52, E2 =1.03 keV, ti=0.19 and 
iV,,=0.38xl022 cm-2. 

3.3. Fitting procedure 

Since the sample is only composed of Seyfert 1 galaxies, 
we expect, in the unification model framework, relatively 
small inclination angles. In the case of NGC 4151, there 
are c laims of a large inclination angle (~60°, Evans et al. 
1993). However this value is rather uncertain. To reduce 
the length of the fitting procedure, we thus decided to fix, 
for all the objects of our sample, the inclination to 30° 
(Nandra et al. |1997D. 



The geometrical normalization of the reflection compo- 
nent R was let free to vary. We recall that a value R = 1 
corresponds to a covering factor of the cold matter to the 
X-ray source of fi = 27r whatever the degree of anisotropy 
of the X-ray emission. 

We added a gaussian to reproduce the iron line complex 
near 6.4 kev. The central energy and width of the gaussian 
were let free to vary. An Iron edge in the range 7-9 keV 
was also added if required by the data, that is for NGC 
4151, NGC 3783. 

We allowed the relative MECS to PDS normalization 
to vary by 5% around the value of 0.86 to acco unt fo r 
the estimated systematic uncertainty (Fiore et al. 1999|) . 
Analogously, we let the LEGS to MEGS normalization free 
to vary o ver the range of acceptable values 0.7-1 (Fiore et 
al. 



1999) 



Finally, in the case of PEXRAV, we fixed the parame- 
ters of the soft excess and/or WA components (i.e. edges 



line, soft power law) to the values presented in section 3.2 



For fits with AG2, we also fixed, in a first step, the column 
density to the best fit values obtained with PEXRAV. A 
second series of fits were obtained with a column density 
free to vary. We checked that there were no significant 
differences between the two sets of results (we generally 
obtain slightly larger column densities and changes of the 
soft temperature kTtb when is let free to vary). We 
will thus only consider the best fit results with free column 
density in the following. 

For computing errors on the parameters of the primary 



continuum (that is kTe, r, kTbi, and R for AG2, and Ec, 
r and R for PEXRAV), we let the column density free to 
vary but we fixed the LEGS/MEGS and MEGS/PDS nor- 
malization, and the Iron line and/or edges to their best fit 
values. Throughout the paper, errors on single parameter 
are quoted at a confidence level of 90 % (i.e. Ax^=2.7) 
unless otherwise specified. 



4.1. 



4. RESULTS 
The PEXRAV model 



The best fit parameter values obtained with the fitting 
procedure described above using the PEXRAV model are 
reported in Table ||. We checked, as far as possible, our 
PEXRAV results with those already published on these 
data. For IG 4329A, Mkn 509, NGG 4151, NGG 3783 
and NGG 5548 our PEXRAV fits are in go od agreement 
with t hose obtained by Perola et al. (1999), Perola et al . 
( |200C| ), Piro et al. (|2000a| E000b| ), De Rosa et al. ( pOOO| ) 
and Nicastro et al. ( |200C ) respectively. In the case of 
NGG 5548, we also checked that a better modelization of 
the W A (with the GLOUDY code used by Nicastro et al. 
2000) does not affect the results on the continuum. We 
note also that, in the case of NGG 4151, we found smaller 
errors for the spec tral ind ex F (AF ~ ±0.01) in compari- 
son to Piro et al. ( |2000a| , where they found AF ~ ±0.1). 
This is simply due to the fact that the latters have com- 
puted the errors using the MEGS and PDS data (above 
3.5 keV) whereas we have used the LEGS, MEGS and PDS 
ones together (thus adding the soft excess parameters in 
the fitting procedure). 

For the PEXRAV model the parameters of the primary 
continuum determined by the fits are F and Ec- The ac- 
tual physical parameters of the Gomptonizing region (also 
reported in Table ||) have been obtained from the spectral 
parameters in the following way. The temperature kT^ is 
simply estimated as kT^ = Ec/2, keeping in mind that 
such approximation roughly holds for t<1. For t ^ 1, 
kTe = Ec/S would be more appropriate. Thus the re- 
ported values can be considered as upper limits to the 
temperature. Knowing the temperature, the optical depth 
can be computed from the spectral index derived from the 
PEXRAV fit usi ng th e following relation (Sh apiro , Light- 
man & Eardley |1976t Sunyaev & Titarchuk |1980| ; Light- 
man & Zdziarski 1987): 



kTeT{l+T/3) 



1/2 



(1) 



This equation is valid for r > 1, and we have checked a 
posteriori that such condition is roughly matched in all 
cases. 



4.2. The slab model 

The fits with the AG2 yield directly and consistently 
the temperature kT^ and optical depth r, which was our 
first motivation for starting this approach. The best fit 
parameters are reported in Table H. It is worth noting 
that another parameter enters the fitting procedure in an 
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important way, that is the temperature of the soft pho- 
tons, fcTfeb- Its values (see Table ||) range, in most cases, 
between 10 and 40 eV and are relatively well constrained 
with errors of the order of 5-50% whereas we have only 
data above 100 eV. For large values of kTbb (i-e. kTtb >30 
eV), the high energy tail of the black body shape may be 
detected in the low energy part of the BeppoSAX data and 
kTbii can then be roughly constrained. We recall that the 
intensity of the seed photon component with respect to 
the comptonized one is also fixed by the fitted parameters 
(mainly r). When no tail is observed, the soft-X-ray data 
impose, at least, an upper limit of fcTbf,- For anisotropic 
geometries, the high energy part of the Comptonization 
spectrum also depends on kT^b, through the location of 
the anisotropy break. The fitting procedure has thus to 
adjust kTbb, together with t, kTe and R, to reproduce the 
soft and hard X-ray data resulting in a complex interde- 
pendence in the fitting procedure, between the soft and 
hard part of the spectra. 

It is worth noting that, in the case of NGC 4151 96L 
and 96II, the best fit values of kTbb are relatively large 
(in fact we have only lower limit, the AC2 working only 
for kTbb < 100 eV). It is true however that in the case of 
the high state 96H, we don't have any LEGS data. Conse- 
quently, as explained above, kTbb is not constrained by the 
soft part of the spectra and large values may be reached. 
But it is surprising that in the low state 96L, we still find 
large kTbb, at least 3 times larger than in 1999. A pos- 
sible explanation is that, in the observation of 1996, the 
presence of a larger reflection component (also found with 
the PEXRAV fits) may hide the anisotropy break, thus 
diminishing the constraints on kTbb, and allowing larger 
values. 

4.3. Comparison between the slab model and PEXRAV 

A purely statistical comparison between the two mod- 
els is not possible since both give acceptable values of x^. 
The slab model has one more parameter and gives gener- 
ally lower values of x^, except for NGC 4151. However, we 
have checked, for this source, that the best fit probabil- 
ity are similar for the two models. Besides, the two best 
fits have residuals without significant features and differ- 
ences between each other. The two models give thus an 
equivalent representation of the data. 

4.3.1. The temperature and optical depth of the corona 

For all objects, the estimated corona temperatures (re- 
spectively optical depths) from PEXRAV are substantially 
smaller (respectively larger) than those inferred with an 
anisotropic Comptonization model. This trend is similar 
to the the result obtained for NGC 5548 by POO. In Fig. ga 
we have plotted the corona temperatures obtained with 
AC2 versus the temperatures obtained with PEXRAV. 
Large differences (up to a factor 8) are found between the 
two estimates. Note however that for IC 4329A, for which 
the high energy data are comparatively better than for 
the rest of the sample, the two models give results which 
are quite close to each other. 

The reason of these differences is relatively simple. 
Within PEXRAV type models, the slope of the power law, 
determined with small errors by the LEGS and MEGS 



data, cannot change, by hypothesis, at higher energies. A 
cut-off around 100 keV is then required to fit the PDS 
data. In Comptonization models, the LEGS and MEGS 
data determine the slope below the anisotropy break. 
Above this break the intrinsic spectrum is steeper and 
thus can fit the PDS data without an additional steepen- 
ing beyond 100 keV, allowing for a larger temperature and 
(consequently) a smaller optical depth to keep, roughly, 
the same power law slope. 



4.3.2. The reflection component 

The reflection normalizations obtained with AG2 clus- 
ter around 1, except for two extreme states of NGC 4151. 
They are, in all cases but one (ESO 141-G55), larger than 
those found using the simple cut-off power-law model. We 
have reported in Fig. |l|b the amplitude of the reflection 
component derived from both models. In some cases we 
obtain differences of factors 4-5. Note that with the ex- 
ception of NGC 4151 for all other sources the slab model 
yields values of R consistent with 1. 

The case of NGC 4151 is peculiar, showing a large vari- 
ation in R between the 96 and 99 observations and an 
indication of a variation of R also between the low and 
high states of the July 96 observation. These changes 
are present irrespective of t he model used in the spectral 
fitting (see also Piro et al. |20004 |2000b| ). However the 
slab model yields higher absolute values, allowing for some 
reflection also in the 99 observation for which PEXRAV 
yields R = 0.01 ±0.1. According to the slab model, the 
corona temperature is very high in the 99 observation, 
which may offer a physical con nection with the low value 



of R as dicussed in section 5.1 



We can also remark that the relatively small reflection 
component and the hardness of the spectrum of this object 
may be difhcult to reconcile with simple corona geome- 
tries (like a plane or a patchy corona for example). More 
compl ex on es (as dynamics coronae for example, Malzac 
et al. |2000| ) may be needed. Another possibility may be 
that the reflecting material in NGC 4151 is highly ionized. 
In this case, fitting with models which do not take into 
account the complex ionization pattern of the reflector can 
severely underestimate the re flection normalizati on (B al- 
lantyne, Ross & Fabian 2000| ; Done & Nayakshin 2001; . 



As an illustration of these results, we have plotted in Fig. 
|2^ and ^ the unfolded best fit spectra of IC4329A and 
NGC 4151 (observation of 1999) obtained with PEXRAV 
(dashed lines) and AG2 (solid lines). We have also plotted 
in each case the reflection hump component to better es- 
timate its contribution to the total spectrum. In the case 
of IC 4329A, the two models give roughly the same result 
at high energy. For NGC 4151 however, large differences 
are expected between PEXRAV and AC2, with a factor of 
2-3 at 300 keV and larger than 10 above 500 keV. 

5. DISCUSSION 

A Comptonizing corona above a passive optically thick 
layer (approximating an accretion disk with negligible in- 
ternal heating) is a complex system, in which the soft 
seed photons for the Comptonization process are pro- 
duced by reprocessing of the Comptonized emission in 
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the underlying passive layer. Such a corona is "self- 
cooled" and, for a given geometrical configuration of the 
system, can only be in equilibrium if the temperature 
and optic al de pth satisfy a definite relation (Haardt & 
Maraschi 199l| ). These relations which essentially cor- 
respond to roughly constant Compton parameters y ~ 

t(1 + r), with y ~ 0.5 and 



kT, 



1 



2 for the slab and hemisphere respec tively , have been ac- 
curately comput ed by Stern et al. (1995; see also for a 
review Svensson |1996| ) . 

In Fig. ^ and ^ the values of t vs. kT^, obtained 
for the 6 sources considered here using the spectral mod- 
els AC2 and PEXRAV, are compared with the theoretical 
relations expected for a Comptonizing region in energy 
balance for a plane (filled rectangles, solid line) and hemi- 
spherical (filled hemispheres, dashed line) configuration re- 
spectively. 

The best fit results obtained with AC2 (Fig. ^) are 
relatively close to the theoretical expectations for the slab 
case even if they tend to fall preferentially above the solid 
line. The data therefore indicate a Comptonization pa- 
rameter larger than for a pure slab geometry, that is a 
more "photon-starved" configuration. A similar result 
was already derived by POO for NGC 5548. For the lat- 
ter source, best fits with a hemispherical Comptonisation 
model were also perfomed. The derived coronal parame- 
ters were in fact in good agreement with the theoretical 
predictions for the hemisphere. However an implausibly 
large value of the refiection component (of the order of 
2) was required, suggesting that the real configuration is 
more complex than these two ideal cases. For the present 
sample we did not perform fits with hemispherical models. 

The values of kT^ and t obtained with PEXRAV (cf. 
Fig. ||d) show the same trend as found with AC2 that 
is larger values of r for smaller values of kT^. They also 
fall above the theoretical expectation for a slab and even 
above the theoretical expectati on for a hemisphere (the 
theoretical curves of Stern et al. 1995 have been extended 



to large r assuming a constant value of the Compton pa- 
rameter for each curve i.e. y ~0.5 and 2 for the slab and 
hemispherical configuration respectively). Therefore this 
set of parameters suggests, for each source of our sample, 
a configuration more "photon-starved" than a hemisphere. 

An important difference with respect to the slab model 
is that the optical depths derived from PEXRAV fits are 
generally larger than 1. The corona should then be opti- 
cally thick reducing or cancelling the effects of anisotropy. 
Therefore also the PEXRAV model has an internal con- 
sistency. However a corona with large optical depth may 
wash out discrete features from the und erlyin g disk (e.g. 



lines and reflection itself, Petrucci et al. 2001) more than 
would be desirable. This problem could be alleviated if 
the corona was "patchy" (and possibly dynamic cf. 
tion[4.3.2D. 



sec- 



We conclude that both a hot, optically thin corona 
with significant anisotropic effects and a less hot, opti- 
cally thick, patchy corona with negligible anisotropy are 
consistent with the available data. 



5.1. Correlations between physical parameters 

A correlation between the reflection normalization R 
and the photon index F has been claimed by Z99 from 
the study, using PEXRAV to model the primary contin- 
uum, of a large number of GINGA observations of Seyfert 
and galactic black hole objects. Yet, the validity of this 
correlation is under debate since R and F are strongly 
corre lated in the fitting procedure (Vaughan & Edelson 
2000). It seems however more significant and less dispute 
in the case of the galactic black hole objects where the 
errors on both parameters are smaller. 

We have plotted in Fig. ^ the reflection normalization 
R versus the photon index F we obtained with PEXRAV. 
In our case, we do not find a clear correlation between 
the two parameters. A Spearman rank-order correlation 
test gives inconclusive results with Ts ~0.3. However our 
sample is biased in favor of objects with hard spectra, to 
ensure a good detection in the PDS instrument. When 
the 13 objects actually observed by BeppoSAX are taken 
into account a stronger correlation is observed (MOO), in 
agreement with Z99. 

In the case of the AC2 model we have no simple way of 
characterizing the spectral shape. Also results are avail- 
able for only 6 sources and we cannot include the other 
Seyfert Is observed by BeppoSAX. Therefore a direct com- 
parison with the Z99 correlation is not possible. Below, 
however, we examine the correlation of R with the corona 
temperature which can be computed for the two models. 

5.1.1. Relation between R or T and the temperature of 

the corona 

Z99 interpret the correlation between R and F in the 
framework of thermal reprocessing models where R is di- 
rectly proportional to the solid angle subtended by the 
cold matter surrounding the corona. The latter solid an- 
gle is not necessarily proportional to the soft photon flux 
reentering the corona, but it is likely that the larger is R, 
the larger is the ratio of soft luminosity to hard luminos- 
ity in the corona, i.e. the smaller is the Comptonization 
parameter y, resulting in softer spectra (larger F). 

It is interesting that the BeppoSAX data can provide the 
temperatures associated with sources with different values 
of R or F, which were not available from the GINGA data 
used by Z99. The result is shown in Fig. ^ where we 
plot R vs. the corona temperature Ec/2 obtained with 
PEXRAV. Here again, in order to have a larger sample, 
we added the best fit values obtained by MOO for the other 
Seyfert Is observed with BeppoSAX. This plot suggests a 
positive correlation between the two parameters, that is, 
the temperature of the corona is larger for larger values 
of the refiection normalization. In fact a plot of Ec/2 
vs. F shows that steeper spectra correspond to higher 
temperatures and (as a corollary) to lower optical depths 
(cf. Fig. |a and|b). 

The slab model analysis yields different trends. A plot 
of R vs. kTe for the latter model shows that sources 
with larger R tend to have lower kT^ (cf. Fig. M), a be- 
haviour opposite to what is found with PEXRAV. The 
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suggested correlation depends however strongly on the 
different states of NGC 4151. 

To summarize, an analysis in terms of PEXRAV indi- 
cates the following correlations: larger reflection compo- 
nent (steeper F) - larger temperature (cf. Fig. ^ while 
the analysis in terms of anisotropic Comptonization yields: 
Zan7er reflection - /oiwer temperature (cf. Fig. 0). 

The latter behaviour is naturally expected in pair free 
Compton cooled coronae (see Fig ^) since the transition 
to lower Compton parameter (i.e. higher Ig/lh as suggested 
by a larger R) implies a decrease in the temperature for 
constant optical depth, as (relatively well) verified by our 
sources. 

The behaviour indicated by PEXRAV could instead 
be understood better in the case of a pair dominated 
corona in pair equilibrium. In such a case, for constant 
hard compactness, the larger the ratio Is/lh: the larger 
the temperature of the corona. This behavior is simply 
imposed by the pair equilibrium. Indeed, an increase of 
the cooling corresponds to a decrease of the number of 
particles in the hard tail of the thermal particle distribu- 
tion. To reach a ne w equi librium the temp erature must 

Ghisellini & 



mcrease. 



(Zdziarski |1985| ; Svensson [1982 

Haardt |1994| ; Coppi [l99g| ) 



A possible limitation of the "pair dominated" interpre- 
tation comes from the high values of the compactness 
Ih required by the physical parameters determined with 
PEXRAV. We have reported on Fig. || the r-Tg relations 
expected, for a pair dominated corona, for different values 
of Z/i, but varying the ratio Ih/h- We have also included 
the observations of MOO. The data require values of It 
in the range 100-1000. Given the (unabsorbed) luminos- 
ity of the sources of our sample (Lo. 1-200 keV of the order 
of lO''* erg.cm^^.s^^), this implies black hole mass upper 
limits of the order of 10^ solar masses (for Z;i=100). This is 
relatively smaller than (but roughly consistent with) black 
hole mass estimations obtained from reverberation map- 
ping or resolved kinematics methods for mos t of the ob- 
jects of our sample (see Gebhardt et al. 



200c and Nelson 



2000 and references therein). The worse case is that of 
NGC 4151, where the PEXRAV best fit values of kT^ and 
r require very large hard compactness (larger than 1000), 
implying a black hole mass of ~ 10^ Mq (reverberation 
mapping met hods e stimate a black ho le ma ss of ^ 10 ''Mq, 
Clavel et al. |l987| ; Ulrich & Home |1996| ). In this case, 
however, the complexity of the soft part of the spectrum 
may lead to an uncorrect modelisation of the continuum. 
We therefore conclude that the pair dominated coronae 
interpretation is also still valid for this object. Other ob- 
servations (with better constraints on the compactness and 
the physical parameter of the corona) are clearly needed 
to confirm or disprove these results. 



6. CONCLUSION 

The aim of this paper was to test anisotropic Comp- 
tonization models over the high signal to noise BeppoSAX 
observations of a sample of six Seyfert 1 (seven observa- 
tions). We use two types of model: a detailed Comp- 



tonization code in slab geometry, which treats carefully 
the anisotropy effects in Compton processes, and a simple 
cut-off power law plus reflection model (PEXRAV model 
of XSPEC). The latter is generally used as a zero order 
approximation to Comptonization spectra. If this is a rel- 
atively good approximation for isotropic geometry, it fails 
to reproduce the real shape of Comptonization spectra in 
anisotropic ones (like a plane, hemispherical or spherical 
hot region above a flat disk). In this case, due to the defl- 
ciency of the first order scattering component towards the 
observer, it better resembles a broken power law with con- 
vex curvature. This leads to strong differences between the 
best fit parameter values obtained with these two models. 
The main results of this work are the following: 

• The data are well fitted by both models and there 
is no statistical evidence for a model to be better 
than the other. Both models give results in agree- 
ment with a X-ray source geometry more "photon 
starved" than the slab case. 

• Our best fits with Comptonization models in slab 
geometry give a temperature generally much larger 
and an optical depth much smaller than derived from 
the power law -I- cut-off fits, using standard Comp- 
tonization formulae. The estimate of the refiection 
normalization is also larger with the slab corona 
model. 

• The two models also lead to different relation- 
ships between physical parameters. For instance, 
PEXRAV indicates a correlation between the re- 
flection normalization and the corona temperature 
whereas the slab corona model suggests an anticor- 
relation. This has major consequences for the phys- 
ical interpretation of the data. Thus, the PEXRAV 
results suggest that the hot corona may be pair dom- 
inated whereas the slab corona model is in better 
agreement with a low pair density solution. 

These results should have observational consequences. 
First, the large differences in the fitted temperatures lead 
to widely different predictions as to the fluxes emitted at 
higher energies, above the PDS range. This underlines 
the need of better soft 7-ray observations, as those possi- 
bly provided by INTEGRAL, in order to directly conflrm 
or disprove the temperatures inferred with anisotropic 
Comptonization models. 

On the other hand, the large values of r predicted by 
PEXRAV would spread out the reflection hump in the 
entire X-ray range, avoiding any detection of this compo- 
nent. The fact that we obtain not negligible values of R 
with PEXRAV may however be reconciled with the large 
values of t if the corona is "patchy" . Indeed, in this case, 
reflection could come from the uncovered part of the cold, 
reflecting material. Part of the reflection could also come 
from the torus. The variability of the reflection component 
may allow to differentiate these two possibilities. Large 
values of r would also modify the profile of the Fe line 
component produced in the disk, leading to an attenua- 
tion and broadening of the line. Forthcoming observations 
with CHANDRA and XMM- Newton are expected to bring 
substantial progress on this issue. 
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Source 


Observation 


Exposure 


PDS 


References 


Name 


Date 


Time (ks) 


Detection 





NGC 5548 
IC4329A 
NGC 4151 



ESO 141-G55 
Mkn 509 
NGC 3783 



Aug. 97 
Jan. 98 
Jan. 99 
Jul. 96L 
Jul. 96H 
Nov. 99 
May/Oct 98 
Jun. 98 



300 
100 
100 
56 
15 
100 
100 
300 



55 
68 
182 
130 
90 
15 
29 
55 



Nicastro at al 
Perola at aL 
Piro et id 
Piro et al 
Piro et al 



2000a, 



2000 



199E 
200M 



2000a, 



2000b 



2000b 



Perola et al. |200d 
De Rosa et al. POOOl 



Table 1 

Sources names, observations dates and exposures of our sample. We have also reported the statistical 

SIGNIFICANCE IN (J OF PDS DETECTION IN 13-300 KEV BAND. We HAVE ALSO REPORTED, AS FAR AS POSSIBLE, THE 
REFERENCES PRESENTING THE FIRST DATA ANALYSIS OF EACH OBSERVATION. 



A^H(1020cm- 




Er ^ 






E2^ T2 






0.55t{J;f 


n 1 +u.ub 

'-'■^-0.04 




^■'-'-'^-0.03 "-^"-0.04 



Table 2 

Best fit parameter of the soft X-ray complex of NGC 3783, i.e. two edges and an emission line (see text for 
details), obtained with PEXRAV. We have also reported the best fit column density, in excess of the galactic 

ONE. 
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Table 3 

Best fit parameters using the PEXRAV model: the photon index F, the temperature fcTe (in keV) and optical 

DEPTH r OF THE CORONA, THE REFLECTION NORMALIZATION R AND THE BEST FIT TOTAL (i.E. INTRINSIC+GALACTIC, IN 
10^^ Cm~^) COLUMN DENSITY FO R EAC H OBJECT. In THE CASE OF NGC 4151, THE COMPLEX NEUTRA L AB SORBER, IS DETAILED 
IN SECTION 3.2.1 . The computation of kTe AND T ARE EXPLAINED IN SECTION 
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Table 4 

Best fit parameters using the Comptonization model of H94: the temperature and optical depth of the corona 
feTe (in keV) and t, the soft photon temperature kTtb (in eV), the reflexion normalization R and the best fit 

TOTAL (i.e. INTRINSIC+GALACTIC, IN 10^^ Cm'^) COLUMN DENSITY FOR EACH OB.IEC T. In THE CASE OF NGC 4151, THE 

COMPLEX NEUTRAL ABSORBER IS DETAILED IN SECTION ^3.2.1 . 
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Figure Captions 



Fig. 1. — (a) Best fit values of the corona temperatures obtained with the cut-off power law + reflection model PEXRAV versus the 
best fit values obtained with the Comptonization model in slab geometry H94. (b) Same as (a) but with the best fit values of the refiection 
normalization. The filled squares represent the different observations of NGC 4151. 

Fig. 2. — (a) Unfolded best fit spectra of IC 4329A with PEXRAV (dashed lines) and H94 (solid lines). For each model, we have also 
plotted the continuum without the reflection contribution and the reflection hump, (b) Same as (a) but for the observation of 1999 of NGC 
4151. 

Fig. 3. — Optical depth r versus temperature kTe from (a) H94 and (b) PEXRAV flts. The theoretical relations between r and kTe for a 
plane and hemispherical Comptonizing region in energy halajjice, for kTijij=5 eV, are shown for comparison in solid and dashed line respectively 
(filled rectangles and hemispheres are from Stern et a|^ |l995| . These theoretical relations have been extended to large optical depth assuming 
the constancy of the Compton parameter, cf. section Bl). The filled squares represent the different observations of NGC 4151. 

Fig. 4. — Refiection R versus photon index F from PEXRAV. The open symbols are from MOO while the solid ones are the results of this 
work. The filled squares represent the different observations of NGC 4151. 

Fig. 5. — Reflection component R versus corona temperature Ec/2 from PEXRAV. The open symbols are from MOO while the solid ones 
are the results of this work. The flUed squares represent the different observations of NGC 4151. 

Fig. 6. — (a) PEXRAV corona temperature Ec/2 and (b) PEXRAV corona optical depth r versus photon index F. The open symbols are 
from MOO while the solid ones are the results of this work. The flUed squares represent the different observations of NGC 4151. 

Fig. 7. — Reflection normalization R versus corona temperature kTe obtained with the slab model. The filled squares represent the different 
observations of NGC 4151. 

Fig. 8. — Same as Fig. ^ but with extra T-kTe curves expected in the case of a pair dominated spherical corona for constant values of the 
hard compactness (in dot-dashed lines. From left to right, Z(,=1000, 100 and 10) but varying the starvness ratio Is/lh- We have added the 
data from MOO in open symbols while the solid ones are the results of this work. The filled squares represent the different observations of 
NGC 4151. 
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